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C3 , Abstract 



00 ' We evaluate three meson doorway mechanisms for nucleon-antinucleon 

annihilation at rest for the first time. Detailed results are presented for the 
final state (jm originating from the 3 Si initial state and for the (ftp channel 
' originating from 1 Sq. The results presented also include the improved 

contributions from two meson doorway states and from the tree diagrams. 
For all the channels considered a consistent explanation of large and small 
f\| ' OZI violations emerges. 



(N 



^ . Recent and accurate data for nucleon-antinucleon annihilation at rest from ex- 

<D \ periments performed at LEAR JIH8[] have challenged our understanding of the 



1 Introduction 



underlying annihilation mechanisms and of the production of mesons with hid- 
den strangeness in particular. Large violations of the OZI rule for special channels 
have been observed. The biggest deviation from the OZI prediction for hadronic 
channels occurs for the (f)ir° final state and has led to speculations about the 
internal structure of the nucleon suggesting a large ss component in the wave 
function [p|— pTlfl . Earlier analysis [12-21] has shown that two meson doorway 



contributions have the correct magnitude to explain the experimental branching 
ratio for this reaction. The present paper extends the preceding calculations by 
including three meson doorway states. Sizable OZI-rule avoiding contributions 
are expected from such intermediate states since the first step, the annihilation 
into three non-strange mesons, represents about one third of the total annihilation 
cross section. 



Based on the results of |22| we shall present a comprehensive effort of cal- 
culating all relevant diagrams involving (non-strange) three meson intermediate 
states leading to the two meson final states 07r° and The corresponding 

two-loop amplitudes have been evaluated with full spin. We have also completed 
the evaluation of one-loop amplitudes (two-meson-doorways) where needed. For 
completeness we report some of the results on other two-meson final states. We 
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(E b , k b ) 



Figure 1: (a) The generic three- meson-doorway diagram; (b) unitarity approxi- 
mation. The notation p n = (E n , k n ) is used for the four-momenta of the particles 
n — a,b, c, x, y, z, v, w in the CMS. 

shall show that the results consistently explain the size of large and small OZI 
violations. 

The paper is organized as follows. In Sec. ^| we describe the three meson door- 
way formalism. A generic case is reported in some detail while technicalities are 
relegated to the Appendices. Section || presents the three-meson-doorway results. 
Updated calculations for the two-meson doorway amplitudes are included and the 
full calculation is compared to the experimental branching ratios. Section f| gives 
the conclusions. 



2 Three-meson doorway mechanisms 

The generic three-meson-doorway diagram is shown in Fig. [I] which also defines 
the notation in terms of four and three vectors. The crosses in Fig. [l|(b) denote 
the unitarity approximation for which all the three s-channel particles are on their 
mass shell. The available data for the annihilation into three mesons relevant for 
the first step of the three-meson-doorway mechanisms are summarized in Table |7| 
of Appendix [A. 2. 



T 



The basic expression for the amplitude without spin is 

i 2 ff g a 9b9c9d d 4 p x d A p y 



(27r) 8 J J (pi — m? x + ie) (p 2 — m 2 + ie) [p 2 — m 2 z + ie) (p 2 v — ml + ie) (p 2 , — m 2 w + ie) 

(1) 

where g a , gb, g c , and g^ are the coupling constants corresponding to the vertices 
a,b,c,d in Fig. |l|(a). In the on-shell or unitarity approximation Fig. |l](b), the 
amplitude T V a has the form 

rp = Wa9b9c9d f d$ 3 (p a ,p x ,p y ,p z ) 
UA " 2(2vr) 5 J (pl-ml + ie)(p 2 w - m 2 w + ie) U 

where <i$ 3 is the phase space of the intermediate three-meson state 

, ^ / \ r/ n C^kr d 3 \t y d?\i z 

d<S> 3 {p a ,P X ,Py,P Z ) = 5{p a -p x -Py-p z )———= r —— (3) 

Lh, x zh/y ihi z 
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The extension of Eqs.f]!]^) to the case of particles with spin is straightforward. 
The details for the vertex spin structure and coupling constants used are given in 
Appendix [A|. The spin formalism for the overall amplitude of Fig.[I] is summarized 
in Appendix [TJ. 

The leading two-loop mechanisms for the reactions pp — > <f)7t° correspond to 
the diagrams in Table [TJ. The <f)ir channel has I G = 1 + and originates from the 
pp( 3 Si) state J PC = 1 . Here we expect that the pirn doorway mechanism is 
important because the pnir system has the largest phase space among all three 
meson states with positive G-parity and is strongly produced in pp annihilation 
(see Appendix [A|, Table 0). The two- loop diagrams in Table [j] all proceed through 
the same doorway mesons (pirn) but differ by the mesons exchanged in the t- 
channel. The tttt subsystem in the intermediate state has total isospin 1 = 
and total angular momentum J = 0. Thus there is no danger of double counting 
when this three-meson doorway mechanism is added to the two-meson pp doorway 
mechanism which is known to be very important [ITS . 



Spin effects have been calculated by introducing scalar invariant functions as 
described in Appendix [TJ. As a cross check, helicity amplitudes for the full ampli- 
tude T have been evaluated, squared and summed. Several charge configurations 
for the intermediate states, see Table [l], add coherently leading to an enhancement 
of the two-loop contributions. The corresponding isospin factors are collected in 



Table |8| of Appendix |B.1| . Because the t-channel particles can reach the mass 
shell, the unitarity approximation acquires a real part. A similar situation has 
been encountered already in the one loop calculation for the two-meson doorway 



mechanism 19 



For the kinematical situations where the t-channel particles are off-shell we 
have introduced monopole form factors 

XI 'Pi X 2 W - {p b +Py) 2 ' 

The parameters have been varied in the range A^, A„, = (1.0 — 1.5) GeV, similar to 
the one loop calculation [O. The form factors reduce the unitarity amplitude by 



about a factor two. For a calculation beyond the unitarity approximation form 
factors for the s-channel doorway mesons must be introduced as well. The corre- 
sponding off-shell contributions are expected to be comparable to the unitarity 



amplitude, similar to the detailed evaluations done in the one-loop case JT9| . 

Turning to the (ftp channel which has I G = 1~ and originates from the pp{ l So) 
state J PC = h , we expect that the tch-it and ujtttt doorway mechanisms are 
important because these intermediate states have the largest phase space among 
all three meson states with negative G-parity. Several states for the nmr system 
are possible in this case which can be classified by the symmetry of the isospin 
wave function. The completely symmetric isospin wave function corresponding 
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Experiment 7.0 ±0.4 



Table 1: Two-loop and one-loop diagrams for the reaction pp^So) — > 07r. The 
corresponding amplitudes T are normalized to the tree-level amplitude g^^ 
from OJ0 mixing. 
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mechanism 




Experiment 4.9 ±0.8 § 



Table 2: Two-loop and one-loop diagrams for the reaction pp^So) — > (pp. The 
corresponding amplitudes T are normalized to the tree-level amplitude g^^p from 
uxf) mixing. 
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to the Young tableau cm has isospin 1 = 1 (another completely symmetric 
state I = 3 is excluded by isospin conservation) J23|, with the space part being 
also completely symmetric. In Appendix |A], the notation is used for the cor- 
responding coupling constant. The state with the mixed symmetry of the isospin 
part EP can also have total isospin 1 = 1, the corresponding coupling constant 
is (see Appendix [A] for details). In this case the 3tt state always contains 
a pion pair with total isospin I n7T = 1 and odd relative angular momentum. It 
is likely that this configuration is saturated by the p ± 7T =F channel (the channel 
p°Ti° is not coupled to the 4>p° because of wrong C-parity). We shall drop such 
three-meson intermediate states to avoid double counting. The completely anti- 
symmetric isospin wave function has total isospin 1 = and cannot occur in the 
annihilation into 07r. Concerning the intermediate state umr, its isospin structure 
is completely determined by the total isospin of the final state 1 = 1 which is 
equal to the isospin of the pion pair. Therefore one can expect that this interme- 
diate state is saturated by the one-loop intermediate state up . For the purpose 
of information the tables show all relevant amplitudes calculated separately. The 
two-loop contributions are fairly sizable. 



3 Results 

The complete results for the (pir and <j)p° channels are summarized in Tables |], 
H, and |3|. The three-meson doorway states shown are the ones leading to the 
biggest contributions. The contributions from other intermediate states have 
been calculated and were found to be negligible p2 |. The 'tree-level' amplitudes 



correspond to the u<p mixing which is proportional to the deviation of the physical 



mixing angle = 37.6° |24| from the the ideal one Gj = 35.3°: 

g^ x = tm(e-e i )-g ar ^ x . (6) 

The on-shell values of the vertex functions occurring in the calculation are 
constrained directly by experimental information. For the annihilation vertices, 
the data are shown in Table |7| and the corresponding coupling constants have been 
parametrized as described in Appendix [A.2| . The remaining vertices have been 
calculated from the measured decay widths of the corresponding mesons and are 



collected in Appendix [A.l|. The values shown in Table H do not include the form 



factors of Eqs. (|4]||) which lead to a reduction by about a factor of two. This should 
be quite a reliable approximation since this reduction is expected to be partially 
compensated by contributions originating from off-shell s-channel propagation, 



similarly to the case of two- meson doorway mechanism ||19| , as mentioned in the 
context of Eqs.fl§@). 

In general, major cancellations between amplitudes corresponding to different 
intermediate states are not likely to occur for the unitarity approximation. Be- 
cause the main contribution in this case comes from the absorptive part of loop 
diagrams, the situation is very different from the well known case of Lipkin can- 



cellations [£5[] where, contrary to our situation, threshold effects are negligible and 
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Reaction / Mechanism 



BR ■ 10 4 



pp — > 07t° (tree-level) 



pp — > IT* if 
pp^ pp^ 
pp —> 



?7T 



)7T 







1-loop) 



97T 







rO 



— > p 7T7T — > 
— > 7Tp 7T — > 07T L 

pp — ► 07T° (1-loop and 2- loop) 
(experiment) 



pp — > (j)TT C 



0.13 

0.5 - 1.4 

0.05-2.0 

0.9-5.1 

0.22 

1.1 

5.9 - 17 

6.5 ±0.7 

7.6 ±0.6 
4.0 ±0.8 



pp 
pp 
pp 
pp 
pp 
pp 



0p° (tree-level) 
(f>p° (one-loop) 



UTTTT 
TXTXTX 



,,0 



9° (1-loop and 2-loop) 
s>° (experiment) 



0.14 
0.1 
0.08 
0.10 

~ 2 

3.4 ±1.0 | 



pp 
pp 
pp 



0?7 (tree-level) 

K*K -> ^r? 

0?7 (experiment) 



0.2 
1.0 

0.78 ±0.21 



Table 3: The branching ratios -Bi? calculated for various doorway mechanisms in 
comparison with the experimental data (in units 10 -4 ). The values marked 1-loop 
and 2-loop contributions correspond to different ways of adding the amplitudes 
coherently. 



specific intermediate states interfere distructively. In the case of the 07r channel 
the situation has been discussed in detail [jnj for the one-loop mechanism. 

For the case of the pp doorway contribution to 07r° we have evaluated the 
full range of the possible coupling constants, see Appendix complementing the 
results in [22 . The corresponding range is indicated in Table [| 



For the OZI-violating final state 07r° the calculated branching ratio in Ta- 
ble H is well within the experimental range. The three-meson doorway con- 
tributions are comparable with the two-meson ones. The range of theoretical 
predictions when adding the amplitudes coherently now easily includes the ex- 
perimental branching ratios while the one-loop results alone are somewhat low. 
For the (ftrj channel, the two- meson-doorway mechanism with the K*K interme- 
diate state has been found to be comparable with the experimental data. The 
two-loop calculations for the (pp final states are also reported in Table |3[ As 
has been mentioned at the end of Section ^ two of the diagrams involve double 
counting with one-loop mechanisms, which is not easily quantified. However, the 
two-loop contributions obviously improve the comparison with the experimental 
branching ratio. 
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At the end of this section we would like to mention that we have also evaluated 
two and three-meson doorway contributions for a number of two meson final states 
(7T7T, pn, pp and pui) without hidden strangeness. It is gratifying to observe that all 
the doorway contributions calculated have turned out to be relatively small when 
compared to the experimental rates for the corresponding annihilation channels. 

4 Conclusions 

We have found that the observed OZI violating enhancement of <j) meson pro- 
duction at rest can be naturally explained by two and three meson doorway 
contributions. In our analysis, there appears to be no need to introduce a large 
ss fraction into the nucleon wave function. The doorway calculations presented 
here are well constrained by experimental information. In the first step of the 
doorway mechanism the annihilation rates into non-strange mesons enter. For 
annihilation at rest these transition rates are well measured. This is particularly 
true for the largest observed OZI violation in pp — > 07r° where detailed infor- 
mation on the spin-isospin dependence of the amplitudes for the annihilation 
pp — > K*K exists. Similarly the meson decay vertices occuring in the second 
step of the one-loop doorway mechanism are directly constrained by the mea- 
sured decay widths. The leading one-loop contribution is thus well determined. 
In the present paper we have shown that the dominant three-meson doorway 
mechanisms (two loops) for pp — > 07r are of similar size as the one loop contribu- 
tions. It is therefore established that the full calculation leaves ample space for 
accommodating any remaining descrepancy with the measured branching ratio. 
At the same time two and three-meson doorway contributions to all the other 
channels involving mesons in the final state are small but not negligible due to 
interference, which again is in agreement with measured branching ratios. Four- 
meson-doorway contributions and higher are expected to be negligible due to 
progressively vanishing probability to rearrange non-strange multi-meson inter- 
mediate states into two-meson final states. We therefore believe that the present 
multiple doorway analysis is qualitatively exhaustive for nucleon-antinucleon an- 
nihilation into mesons at rest. Extending these calculations towards higher 
energies seems desirable. However, the experimental information on the energy 
dependence of the production of the intermediate states is far less detailed and 
large uncertainties in the corresponding predictions appear to be unavoidable. 
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A Coupling constants 



In this appendix, we collect the coupling constants for pp annihilation at rest 
into various two and three-meson annihilation channels and the required coupling 
constants for meson decays. For most cases, the couplings can be expressed in 
terms of partial decay widths which are known from experiment and provide a 
model-independent input to the calculations of the doorway mechanisms. 



A.l Vertices for two-particle decays 

The amplitudes corresponding to the transitions a — > b + c involving pseudoscalar 
fields and vector fields with minimal number of derivatives have the following 
form in momentum space 

(4>b4> c \T\V a ) = g V4>4> e a - (p b - p c ) (Al) 

(V b <f> c \T\V a ) = gvV4>^uapP^ U aPt4 (A2) 

(V b V c \T\V a ) = g ( y ] p a ■e b e a -e c + g ( y ] p a -e c e a -e b + g { y ] p b -e a e b - e c (A3) 

where p a , p b , p c are the corresponding four-momenta, e a , £ b , £ c are the polarization 
vectors and e^^fi j s ^he totally antisymmetric Levi-Civita tensor. The decay 
widths T a ^ b+C are related to the corresponding coupling constants g abs by 

r_ 6+c = g ^SM^i ( A4) 

J (ml - (m b + m c ) 2 ){m 2 a - (m b - m c ) 2 ) 

hc = 2^ (A5) 

where k bc is the CMS momentum of the particles b ans c, m a , m b , and m c are the 
corresponding particle masses and f(k bc ) are the spin- weight functions defined in 
Table |. 

The following coupling constants for the meson decays were used in the present 
calculations: g p7Tn = 6.00, g^xR — 4.6, gK*Kn = 5.54, g^ = 1.86 GeV -1 . The 
numerical values of the coupling constants for the two-meson pp annihilation are 
summarized in Tab. ^| together with the corresponding experimental branching 
ratios. For the sake of convenience, the pp annihilation coupling constants are 
normalized to the partial widths of the ground state of the pp atom. These partial 
widths are related to the corresponding annihilation cross sections a pp ^ b+c by 

T pP -> b+c = (va pP ^ b+c ) v ^ \^is(0)\ 2 (A6) 

where v is the relative velocity and l^is^O)) 2 = g7i . p is the probability density 
for the IS atomic state at zero separation between p and p. This gives the 
following relation between the coupling constants listed in Table ||] and the S- 
wave annihilation amplitudes g p+p ^ b+c at zero energy: 

- a 2 KMo^_ 2 *H (A7) 
y a ^ b +c — y P +p^ b +c m — y P +p^ b + c v-™- 1 > 
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Reaction a —>■ b + c 



Spin-weight functions f(k. 



be ) 



0+ - 






o- - 


-> 1-1- 




1- - 


-0*0* 




0±- 


- i-o± 




1- - 


-> i~o- 




1- - 


-> 1-1- 


(a) 


1- - 


-> 1-1- 


(b) 



3m 2 



m\{ml + 1m\ + 2m 2 )fc 2 c 



3m 2 m 2 



2fefe c / g _|_ "'<,"/., 



mlkl 



Table 4: The spin-weight functions f(kb c ), Eq. (|A4|) , for two particle decays. For the 
three- vector-meson vertex Eq. (|A3|) , the two cases correspond to the situations consid 
ered in Section [C]: (a) gi 12 ^ = g\P = —gv , gv = and (b) gv 
both assuming m\, = m c . 



(2) 

9v 



0, g v 3) + 0, 



Process 


Ref. Branching ratio -Bi? 


Qabc 




§ (3.4 ±1.0) • 10~ 4 


g a ^ P = 8.71 • lO" 4 


pp(Zzg.) -> up^+n- 


|26| (2.26 ± 0.23) • 10" 2 


g a ^ P = 2.77 • io- 3 






7fl (6.5 ±0.6) • 10- 4 
| (5.5 ±0.7) • 10- 4 
1 (7.57 ± 0.62) • 10~ 4 
3 (4.0 ±0.8) • 10- 4 


g a ^n = 3.43 • 10- 4 


pp{S ^UJTT 


1 


| (5.7 ±0.5) • 10~ 3 
Z7J (5.2 ±0.5) • 10~ 3 


= 7.59 • io- 4 


pp^ pp 




12] 2.4 ■ 10~ 2 


g% = 8.21 ■ IO" 4 

g f) = 9.79 • 10~ 4 GeV" 1 


pp -> If AT, if if 


|2^ 0.23 • 10~ 2 


g K » R = 7.0 ■ 10- 4 GeV" 1 



Table 5: The experimental branching ratios for the two- meson pp annihilation at rest 
and the corresponding coupling constants used in the present calculations. The coupling 
constants are normalized to the total width of the atomic (pp)is state T( p p) ls = 1 keV 
(the singlet-to-triplet ratio 1 : 3 is assumed for the pp spin fractions). 
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A. 2 Vertices for three-meson annihilation 



In the case of three-particle transitions a^x + y + zwe consider reactions of 
the following types 

0" -> (T(T(r 
0" -> 1 _ 0~0~ 

l - -> i _ o - o- 

The corresponding amplitudes with a minimal number of derivatives read 

(4>x4> y 4>z \T\ 4> a ) = g PPP p (A8) 
(MyV t \T\<f> a ) = g pppv e^ af3 pftfafieZ (A9) 

(0a;02/K|^|K) = Qvppa 6 *' 6 * ( A10 ) 

The coupling constants g axyz are related to the three-body transition widths 

^- a^xyz by 

Va^xyz = 9 T\r I W rf$ 3 . (All) 

where u> are the kinematical factors given in Table |G] and the three-body phase 
space d$3 is defined by Eq.(|3|). 

For the pp(0~) — > imn vertex we have two cases considered in Sect. 0. The 
vertex that is completely symmetric in the isospin of the 7T7r7r system corresponds 
to Eq . flA8|) with the coupling constant gf % . The vertex with the mixed symmetry 
is given by 

(n + n\-\T\pp(J = 0~,I=l)) = g£ K (p -p + +p -p„-2p + -p + ) (A12) 

where p + , po, p~ are the four- momenta of the corresponding pions. 

The numerical values of the coupling constants for the three-meson pp anni- 
hilation are summarized in Tab. together with the corresponding experimental 
branching ratios. 



B Evaluation of Amplitudes with Spin 

The on-shell approximation Fig. |I](b) for particles with spin leads to the following 
expression for the covariant amplitudes Tjja replacing the spinless case of Eq.(0) 

rp ig a gbg c gd r w(p x ,p y ,p z )d$ 3 (p a} p x ,p y} p z ) 

UA ~ 2(2vr) 5 J (pl- m l + ie )(pl-ml + it) 1 ' 

The evaluation of the spin-weight functions w(p x ,p y ,p z ) can be done on the am- 
plitude level using a decomposition of the spin functions into covariant tensor 
structures built of the external momenta p a and pb- The corresponding computa- 
tions have been done using symbolic codes written in MAPLE and were verified 
by hand for several cases. As a further, independent check we have determined 
helicity amplitudes in an explicit polarization basis, see e.g. Ref. ||32|| , and calcu- 
lated the sum of the squared helicity amplitudes. 
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Reaction 

a — > x + y + z 


Spin-weight factor 


o- -> o- o- o- 
o- -> i-o- (T 

\- -> i-o- o- 


1 

- (k.k y ) 2 ) 

n i K 

1 H - 

3 m 2 . 



Table 6: Spin-weight functions w for different three-body final states. The k x 
and k y are the 3-momenta of particles x and y in the CMS. The interaction terms 
are defined in Eqs.( [A^ - |A10|) . 



Process 


Ref. Branching ratio -Bi? 


1 9 abed | 


pp(liq.) — > c<J7r7r 

-> ^ ) -> 


[2G| 


0.066 ±0.006 
0.0226 ± 0.0023 


Opp^So^un™ = 0-21 GeV -3 


Pp(S) — > TT°7Y + 7Y- 

pp{S -> ^o) -> (7r 7r + 7r") p/l . sp . 
-> ^o) -> 


H 0.066 ±0.008 

0.066- (0.083 ±0.029) 
0.066- (0.014 ±0.006) 


4(^0)^ = 0-015 

flSpsbw^ = °- 0056 


pp(liq.) — ► 7r°7T + 7r _ 


[0 0.069 ± 0.004 




pp(liq.) —>■ 7r°7r°7r 

Pp(S -> -> 7T 7r 7T 


[§ (6.2 ± 1.0) - 10" 3 
6.2 • 10" 3 ■ 0.54 


4(^0)^3. = 0.010 


pp(5 -> 3 5i) -> p^ 7r + 7r -a^ 7r+7r - 




7.61 ■ 10~ 2 ■ 0.50 


9pp(. 3 Si)->p7rx — 0.043 



Table 7: The experimental branching ratios for the three-meson pp annihilation 
at rest and the corresponding coupling constants used in the present calculations. 
The coupling constants are normalized to the total width of the atomic {pp)\s 
state r(pp) ls = 1 keV and the singlet-to-triplet ratio 1 : 3 is assumed for the pp 
spin fractions. 
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B.l Invariant Amplitudes and General Tensor Decompo- 
sition 

The covariant integrals over internal momenta p% and p^ can generally be ex- 
pressed in terms of linearly independent tensors constructed of the external mo- 
menta, p^ and p b multiplied by invariant amplitudes. For the one-loop diagrams 
(see Tables [I] and 0) we introduce the following notation 

(p£> = / f(p*,Pa,p b )p»d% (B2) 

where the scalar function f(p x ,p a ,Pb) is defined by a direct calculation of the one 
loop diagram. The general form of this integral is given by 

(jgHifpS + ^rf (B3) 
where the coefficients l[ and can be found straightforwardly: 

j(i) = / kbEx ~ E bk x z x . . fi4 . 

4 1} = <^> • (B5) 

Here the notation (...) is defined similarly to Eq. (|B2|) with the same scalar func- 
tion f(Px,Pa,Pb)- The second order expressions in the one-loop case contain in 
addition two products of p^ and p b and the metric tensor g^ u : 



where 



(p» = '? } ST + 4 2) pX + 4 2, (p!X + pW) + if PK (B6) 



/« = (-ii^iM) (B7) 

(2) , 2 (fcg£| + E 2 x )k} - 4 E t k t E x k A - (1 - 3 z|)mgg v 

' 2 = ( 2^ > m 

r (2) {{l-3zl)E b k x + 2k b E x z x )k x 

h ~ { 2^ > (B9) 



/(2) = (l-3^)fc 



o,2 - )• ( B10 ) 



4 ~ X 2*g 



In the case of two-loop amplitudes we have more internal degrees of freedom, 
and Eq.(|B6|) is replaced by 



{P X Py)=J f(Px,Py,Pa,Pb)p^Pyd%d%. (Bll) 

The tensor decomposition for this expression takes the following form 

<P» = ^ ST + A XV) P>a + 4 XV) (p£# + PM + I [ : V) P^Pl (B12) 
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where 



fry) 



n 



1 2 



r( x y) 



r( x v) 



. i^Z Z X Zy^k X ky^ 



(B13) 



,2 (k x k y z x Zy + E x E y )k 2 4 E b k b E x k x z x [Z 3 z x Zy)m 2 k x ky 

\ ' o— 91-2 ■ TP i4 J 



2 m;[fc£ 

((Z - 3 z x z y )E b k x + 2 k b E x z y )ky 
2m a k 2 ' 



y Z Z x Zy) k x ky 

2kl 



>■ 



(B15) 
(B16) 



A generalisation to tensors of higher rank is straightforward. In particular, the 
following expressions for a rank three tensor appear in the two-loop diagrams: 



r(yxx) 



n 



(yxy) 



,(z x Z — z y )k 2 x ky ((Z — z x Zy)E x k x k y — (1 — z 2 x )Eyk 2 x )k b /T317 x 

( 2 2E b } (B17) 

(z x - z y Z)k x k 2 {{Z - z x Zy)E x k x k y - (1 - z 2 y )E x k 2 y )k b 
{ 2 2E b 



The different topologies of the two-loop diagrams involving vector and pseu- 
doscalar mesons are shown in Fig.|2]. The corresponding complete spin weight 
functions Eq.(|Bl|) for the diagrams A-D are listed below. 





Figure 2: Different three-meson doorway processes involving pseudoscalar parti- 
cles (thin lines) and vector particles (thick lines). 



yk X j ky, Z Xj Zy ) = (k b E y + E c k y Zy)(k b E x - E b k x z x ) + (B19) 

+ ~ — {— ({Z X + ZyZ)E c ky + (Z X Zy + Z)k b Ey) El + (z x z y + Z)E 2 k b E x \ k x ky+ 



2m a k b 



z x + z y Z)k + 2m z x ) E x k x - i(z y + z x Z)k x + 2m x z y ) E y k 



k 2 - 
K b 
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+ 



(z 2 + l)k 2 + 2m 2 + ( z * z v ~ Z ) k * k y) E * + 
+ Uz 2 x + l)k 2 x + 2m 2 x + (z*zy - Z)k x k y ) E 



K + 

z y + z x Z)E 2 c k 2 x k y + {(z x + z y Z)E x k x k 2 - (z y + z x Z)E y k 2 x k y + 

+ ((2k x Z X kyZy + (Z 2 + l)k 2 y )E x + (2k x Z x kyZy + + l)A£)£ w ) fcj) £ c + 

+ (2(k x z x + k y Zy)E x Ey + ((z 2 + l)k 2 + 2m\ + (z x z y - Z)k y k x )k x z x ^j k 2 h + 
+ {{z x z y + Z)E x E y k x ky + ((z^yZ - z 2 x - z 2 - Z 2 )k 2 - (z 2 x + l)m 2 y )k 2 A k b 



E h 



2(k x z x + k y Zy)E x Ey + ((z 2 x + 1) A£ + 2m* + ( Z x Zy Z^jk X ky)kyZy \ 

- ((z^ + Z)E x E y k x ky + ((z^Z - z 2 x - z 2 y - Z 2 )k 2 x - (z 2 y + l)m 2 x )k 2 y ) k b ] E c ) 



wb ( k x , ky , 2^ , 2^ t 



2(/{ 2) +/ 7 



m a k b 



(Pa -PyPc-Py-Pa- Pc mfj + 



r( x v) i rC^y) b 



, I (Pa • Pj/ Pc-Px-Pa- Pc Px ■ Py) (B20) 

m a Kb/ 

W C (k x , k y , Z x , Z y ) = If' {p c -PwPy-Pv-Pc- Pv Py ■ Pw) + I\ ] Py ' Pw ~ H Py ' Po + 

r(y%y)„ ■ 



+ (l^ yxx) p c ■ p w - I { 7 yxy) p c ■ p, 



(B21) 



(1) 



(rn| - m 2 c )(m 2 y -m 2 c + 2(p a - p x ) ■ (p c - p y ) 
m 2 



~(Pc + Py) 2 + 2(p a - Px) ■ (Pc + Py)) 



(B22) 



The isospin factors resulting from the summation over all intermediate states 
in the diagrams of Tables [TJ and are listed in Table [5]. 



Mechanism 




Isospin factor 


Pp( 3 Si) - 


-> npn - 




1 


Pp{ 3 Si) - 


-> pirn - 


-> 07T° 


i 

3 


PP^Sq) - 


-> TITITI - 


- 0P° (A) 


5 


pp^So) - 


-> TITITI - 


- 0P° (B) 




PP^Sq) - 


-> 7T7ra; - 




PI 
V 3 


pp^Sq) - 


-> 7TC^7T - 




[2 
V 3 



Table 8: The isospin factors corresponding to the three-meson door-way mech- 
anisms pp — > xyz — > (j)X. 
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C Re-evaluation of pp 



PP 



7T 



For the pp diagram the effective Langrangian is not unique since the annihila- 
tion vertex pp — > pp allows for two invariant couplings with minimal number of 
derivatives. As in [|13j, these are denoted by 



T (i) 
1 pp^pp 

T (2) 
1 PP->PP 



91 \^pp ' kpp ■ Ey Epp • Ey kpp • Ex) 

92 £pp ' k x E x ' Ey 



(CI) 
(C2) 



where k x and k y are the four-momenta of particles x and y and E p p, E a 



and E y are 



the polarization vectors of the corresponding particles. In |13| the two cases were 
calculated separately. Here we consider the coherent sum of both amplitudes 



T 



Tpp^pp cos 9 + Tpp_+pp sin 8 



,(2) 



PP^PP 



(C3) 



where the total strength is normalized to BR(pp 
ical estimate in 



12|]. Figure 



pp) = 2.4 % from the theoret- 
shows the result of the calculations of the branching 
ratio for the </>7r° final state (including finite width effects) for the full space of 
parameters g\, #2- To check the self consistency of the door-way calculations 
the contribution of the pp intermediate state to the un° production has been 
evaluated as well. Since the um° channel is not suppressed on the tree-level, it is 
gratifying that the one-loop corrections turn out to be small in comparison with 
the experimental data for this channel. At the same time, the contribution to the 
OZI suppressed channel 07r° is very significant for a broad range of the relative 
strength of gi and g 2 - 

(a) (b) 



m 



A=1.6 GeV 
X=1.4 GeV 




or 
m 




6 (rad) 



9 (rad) 



Figure 3: The branching ratios corresponding to the door-way pp mechanism for (a) 
07T° and (b) luit production from pp( 3 Si) annihilation. The parameter i? defines the 
relative strength of the two couplings in the pp — * pp vertex, see Eqs. CI and C2, 
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